Bloch oscillations in tight-binding lattices with defects 
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Defects in tight-binding lattices generally destroy the onset of Bloch oscillations (BOs) and the 
periodic self-imaging of the wave packet due to the lack of an equally-spaced Wannier-Stark ladder 
spectrum. Here it is shown that localized and extended defects in the lattice can be engineered to 
be transparent for BOs. Such lattices are synthesized from the defect-free lattice by the technique 
of intertwining operators generally employed in supersymmetric quantum mechanics. The energy 
spectrum of the synthesized lattices differs from the Wannier-Stark ladder of the defect-free lattice 
because of the missing of pairs of resonances in the ladder, thus ensuring the persistence of BOs. 
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I. INTRODUCTION 



Bloch oscillations (BOs) represent a fundamental 
coherent transport phenomenon originally predicted 
for quantum particles in periodic potentials driven by 
an external dc force. BOs have been observed in a 
wide variety of quantum and classical physical systems, 
including semiconductor superlattices Jl|, ultracold 
atoms [2[, Bose-Einstein condensates J3|, arrays of 
evanescently-coupled optical waveguides [4|, [5[ .photonic 
superlattices [6], and acoustical superlattices In the 
absence of dephasing and scattering processes and for a 
negligible Zener tunneling, the particle accelerated by 
the external force undergoes an oscillatory (rather than 
a translational) motion owing to Bragg scattering off 
the lattice. The periodicity of the motion is basically 
ascribed to the transition of the energy spectrum from 
a continuous band with delocalized Bloch eigenstates 
in absence of the external force, to a discrete ladder 
spectrum with localized Wannier-Stark (WS) eigenstates 
when the external force is applied and Zener tunneling 
is negligible. A single-band tight-binding lattice in 
one spatial dimension provides the simplest model to 
describe the onset of BOs and the formation of WS 
localized states [H, The inclusion of lattice disorder, 
nonlinear it ies, inhomogeneities or defects has generally 
a detrimental effect on BOs. Several authors have in- 
vestigated the onset of BOs in generalized tight-binging 
lattice models. Among others, we mention BOs in 
nonlinear lattices [lOj and in aperiodic lattices with 
long-range correlated disorder BOs with spatially- 
inhomogeneous dc fields [l2|, BOs in quasicrystals [l3j |, 
BOs for interacting bosons [14], BOs in disordered 
lattices with interparticle interaction [l5j |. and BOs in 
lattices with inhomogeneous intersite couplings [16]. The 
introduction of localized or extend defects in the lattice, 
for either the intersite couplings or the site energies, is 
generally expected to destroy the WS ladder spectrum 
of the defect-free lattice, and thus the onset of BOs. 
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In this work it is shown that localized and extended 
defects in the lattice can be engineered to be transparent 
for BOs. Such lattices are synthesized starting from 
the defect-free lattice Hamiltonian by application of 
the technique of intertwining operators, successfully 
employed in supersymmetric quantum mechanics [13, LL8j 
to add or delete energy states of a given Hamiltonian. 
In our case, the spectrum of the synthesized lattices 
differ from the Wannier-Stark ladder of the defect-free 
lattice because of the missing of pairs of resonances in 
the ladder, which ensures the persistence of BOs and 
self-imaging of a wave packet. The paper is organized as 
follows. In Sec. II, the technique of intertwining operators 
for spectral engineering of tight-binding Hamiltonians is 
presented. In Sec. Ill, the technique is then applied to 
synthesize lattice models with defects that support BOs. 
Finally, Sec. IV outlines the main conclusions. 



II. TIGHT-BINDING LATTICE ENGINEERING 

Let us consider a one-dimensional tight-binding lattice 
described by the Hamiltonian 

H = ^ Kn (\n - l)(n\ + \n)(n - 1|) + ^ V n \n)(n\ (1) 

n n 

where |n) is a Wannier state localized at site n of the 
lattice, n n is the hopping rate between sites \n — 1) and 
|n), and V n is the energy of Wannier state |n) in presence 
of the external applied force. To ensure the Hermiticity 
of the hopping amplitudes n n and site energies V n 
must assume real values. For a defect-free lattice and 
a homogenous force, the intersite coupling n n is inde- 
pendent of n {n n = k) and the site energy V n increases 
linearly with n, i.e. V n = -Fn, where F is the potential 
gradient. The corresponding tight-binding Hamiltonian 
will be denoted in the following by Ho and referred to 
as the WS Hamiltonian. As is well known, its spectrum 
is purely discrete and the allowed energy levels are given 
by Ei = IF, where I = 0,±1,±2, ... [9]. Our goal is 
to synthesize a lattice Hamiltonian, of the form given 
by Eq.(l), with an energy spectrum which differs from 
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that of the WS Hamiltonian Ho because of the missing of 
some resonances in the WS ladder spectrum. This spec- 
tral engineering problem can be solved by the technique 
of intertwining operators generally employed in super- 
symmetric quantum mechanics [17| . which is also valid 
for matrix Hamiltonians (see, for instance, [Hj and ref- 
erences therein). The technique is briefly described in a 
rather general manner in this section. 
Let us indicate by Hi the tight-binding Hamiltonian de- 
fined by Eq.(l) with hopping amplitudes and site energies 
given by {ftn\ Vn 1 ^}, and let us assume that n n — >• k as 
n —> ±00 and that \V n \ is bounded or diverges with an 
algebraic law \V n \ ~ n as n — >> ±00; note that such as- 
sumptions are satisfied for the WS Hamiltonian Ho- Let 
us indicate by /ii, \i2, ^3, ••• the point spectrum of Ho 
and let l^ 1 )) = J2 n $n*\n) be the (proper) eigenfunction 
of Hi corresponding to the energy /ii, i.e. jl9| l 
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with \4&\ —> for n — >• ±00. It can be readily shown 
that, provided that (ffi does not vanish, the following 
factorization for Hi holds 



Hi = QiKi+m 



where 
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Let us then construct the new Hamiltonian H2 obtained 
from Hi by interchanging the operators IZi and Qi, i.e. 
let us set 



^2=^iQi+/ii. 



(10) 



Using Eqs.(4-9), from E.(10) it can be readily shown that 
H2 describes the Hamiltonian of a tight-binding lattice 
[i.e., it is of the form (1)] with hopping amplitudes and 

site energies {^n\Vn 2 ^} given by 
*( 2 ) 
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An interesting property of the new Hamiltonian H2 is 
that its energy spectrum is the same as that of Hi, apart 
from the lack of the discrete energy level E = In 
fact, let us indicate by \^e) — ^2 n ^n(E)\n) a proper (or 
improper) eigenfunction of Hi with energy E. Note that, 
if E belongs to the point spectrum of Hi, \tp n (E)\ — >> 
as n — >• ±00, whereas if E belongs to the continuous 
spectrum of Hi, \ip n (E)\ remains bounded as n —> ±00. 
Let us first assume that E ^ Using the factorization 
(3) for Hi, the eigenvalue equation Hi\i/je) — E\i/je) 
reads explicitly 



QiIIi^e) = (E-fii)\^ E ) 



(13) 



from which it follows that 7Zi\^e) 7^ since E ^ \±i. 
Applying the operator to both sides of Eq.(13), one 
obtains 

niQi$ E ) = (E-fii)$ E ), (14) 

i.e. H2\^e) — E\t\) E ), where we have set \i(je) = T^i\^e) 
or, explicitly [see Eq.(5)] 

ME) = r^ME) + f«^-i(£). (15) 

Therefore, \i^e) is an eigenfunction of H2 corresponding 
to the energy E. Also, from Eqs.(6), (7), (15) and from 
the assumed asymptotic behavior of Kn ^ and it fol- 
lows that \^e) is a proper (improper) eigenfunction of H2 
in the same way as \^e) is a proper (improper) eigenfunc- 
tion of Hi. In a similar way, one can show that any eigen- 
value E of H2, belonging to the continuous or to the point 
spectrum (with E 7^ /ii), is also an eigenvalue of Hi- 
Therefore the continuous and point spectra of Hi and 
H2 do coincide, apart from the E = fii eigenvalue which 
needs a separate analysis. For E = /ii, the (proper) 
eigenfunction of Hi is by construction l^ 1 )) [see Eq.(2)], 
which satisfies the condition ^il^ 1 )) = 0. On the other 
hand, from Eq.(10) it follows that the eigenvalue equa- 
tion Hi\ty) = is satisfied for either = IV^ 1 ^) or 
\^) = |^ (2) ), where Qi|^ (1) ) = and Qi|^ (2) ) = 
The equation QilV^ 1 ^) = reads explicitly 
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Using the expressions of and q n ^ given by Eqs.(6-9), 

the difference equation (16) for can be solved in a 
closed form, yielding 



V4 J 
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In view of the asymptotic behavior of (ffi and tvfi as tl 
±00, it turns out that ipn^ is unbounded as n —> ±00, i.e. 
it is not an eigenfunction (neither proper not improper) 

of H2- Similarly, as (j)^ — » as n — » ±00, the equa- 
tion Qi|V> (2) ) = for |^ (2) ) reduces to Eq.(16) in 



the asymptotic limit n —> ±00, and thus also iph is un- 
bounded at n — » ±00. Therefore, none of the two linearly 
independent solutions IV^ 1 ^) and \^ 2 ^) of the second- 
order difference equation II2W) — ViW) are bounded, 
i.e. ii\ does not belong neither to the point spectrum nor 
to the continuous spectrum of 

The factorization method can be iterated to construct 
new Hamiltonians H4, H5, ... whose energy spectra 
differ from that of Hi owing to the missing of the dis- 
crete energy levels {/ii,/i 2 }, {/ii, /i 2 , M3}, {fJ>i, M2, M3, /m}, 
... For instance, to construct the Hamiltonian let |^) 
be the (proper) eigenfunction of Hi corresponding to the 
energy E = /i 2 , and let us set \<j)^) = Hi\ip). From the 
previous analysis, it follows that is the (proper) 

eigenfunction of H2 corresponding to the energy E = /i 2 , 
i.e. 

4 a) ^i + + = (is) 

Let us then construct the new operators IZ2 and Q 2 , 
defined as in Eqs.(4-9) but with and replaced 

by cjffl and «n\ respectively. The factorization H2 = 
+ then holds. Reversing the order of the 
1Z and Q operators, one obtains the new Hamiltonian 
H3 = 7^2 Q2 + M2, which possesses the same energy spec- 
trum of Ho, except for the missing of the two energy 
levels E = jii and E = /i 2 . The hopping amplitudes 
fti 3 ^ and site energies of the lattice corresponding to 
the Hamiltonian H3 are given by Eqs.(ll) and (12), with 

4>n \ Tn \ and Vn^ replaced by (/4 2) , rffl, and 
(2) 

Vn , respectively. It should be noted that the technique 
of intertwining operators so far described could gener- 
ate lattice Hamiltonians with complex-valued hopping 
rates n n or site energies V n ; for instance, the hopping 
amplitudes of H2 might become complex-valued when 
cj)n^ does not have a defined sign [see Eq.(6) and (11)]. 
This situation, corresponding to a non-Hermitian lattice 
with real- valued energy spectrum, will not be considered 
in this work. It should be nevertheless observed that 
iteration of the intertwining operator technique could fi- 
nally restore the Hermiticity of the lattice Hamiltonian, 
in spite some of the intermediate Hamiltonians are not 
self-adjoint. This is precisely the case of our interest that 
will be discussed in the next section. 



III. BLOCH OSCILLATIONS 

Let us consider the WS Hamiltonian Ho, defined by 
Eq.(l) with K n = k and V n = —Fn. As is well known 
(see, e.g., |8|), Ho has a purely point spectrum (the WS 
ladder spectrum) with energies [see Fig. 1(a)] 

Ei = IF (19) 

and corresponding localized eigenstates 

\u^) = J2Jn + i(l)\n), (20) 

n 
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where 7 = 2n/F and I = 0, ±1, ±2, ... is the quantum 
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FIG. 1: (color online) (a) Schematic of the energy levels (WS 
ladder) of a WS Hamiltonian Ho (n n = k, V n = —Fn); the 
energy level spacing is uniform and equal to F. (b) and (c) Pe- 
riodic breathing and oscillatory modes in a WS Hamiltonian 
corresponding to F = 0.2, ft = 1 and to single-site excitation 
[c n (0) = S Uj o in (b)] and broad Gaussian wave packet excita- 
tion [c n (0) = exp[— (n — 10) 2 /64] in (c)] of the lattice at initial 
time t — 0. Note the periodic self-imaging of time evolution 
at multiplies of the BO period Tb = 2tv/F. (d) and (e): same 
as (b) and (c), but for the WS Hamiltonian with a defect in 
the hopping rate (« n = 1 for n / 1, ki = 1.5). 

number [9[ . Note that the WS state with quantum num- 
ber I is localized at around the lattice site n = —I. Owing 
to the equal spacing of WS modes, the temporal evolu- 
tion of any wave packet is periodic and self-imaging is at- 
tained at times multiples of the Bloch period Tb = 2tt/F. 
This effect is clearly visible by observing the breathing 
or oscillatory modes [8[ corresponding to either an ini- 
tial single-site or broad-site excitation of the lattice, as 
shown in Figs. 1(b) and (c). The introduction of some de- 
fects in the lattice, in either the hopping amplitudes K n or 
site energies V n , generally breaks the self-imaging prop- 
erty of the lattice [see, as an example, the simulations 
shown in Figs. 1(d) and (e)]. In this section we aim to 
construct tight-binding lattices with defects in which the 
self-imaging phenomenon of the WS (defect-free) Hamil- 
tonian Ho is maintained. This goal can be achieve by 
the application of the intertwining operator technique de- 
scribed in the previous section assuming as the starting 
Hamiltonian Hi the WS Hamiltonian Ho, i.e. Hi = Ho. 
The main idea is that any new Hamiltonian H 2 , H3, H^ 
... obtained from Ho by successive application of the 
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FIG. 2: (color online) Behavior of (a) the square of hop- 
ping amplitudes [n n ] , and (b) of site energy offset AVd = 

(2) 

Vn + Fn for the Hamiltonian Hi synthesized from the WS 
Hamiltonian Ho by removal of the WS resonance fi\ — (see 
the inset). Parameter values are k = 1 and F — 0.6. 



intertwining operator technique has a spectrum which 
differs from that of Ho by the missing of some of the 
WS resonances, and thus a periodic temporal dynamics 
of the wave packet is maintained with the same period 
Tb of the original WS Hamiltonian. As a first step, let 
us construct the Hamiltonian H2 by removing from the 
spectrum of the WS Hamiltonian Ho one WS resonance, 
for instance the one with energy \i\ = corresponding 
to the quantum number I = in Eq.(19). According 
to Eqs.(ll) and (12) with k n 1] = k, = -Fn and 

(j)n^ = J n (l), the hopping amplitudes and site energies 
of the lattice associated to H2 read explicitly 



' Jn(l)Jn-2(l) 

4 2 -i(7) 



VP = -Fn + /: + ? (7) -.- 4(7) 



Jn{l) Jn-l(l)' 



(21) 
(22) 



A typical behavior of [^ 2) ] 2 and AV^ Z) = -V r 
Vn + Fn, as predicted by Eqs.(21) and (22), is shown in 



,(2) 



r{2) T/ (l) 



Fig.2. Note that and AVd settles down to con- 

stant values as n —> ±00, so that far from the defect near 
n = the lattice described by H2 behaves like a defect- 
free WS lattice. Note also that, owing to the asymptotic 
behavior of Bessel functions J n at large indices, one has 



,(2) 



(2) 



0asn4 +00 but AV^ 2) 2^/7 = F as 
n — » —00 [see Fig. 2(b)]. Unfortunately, due to the os- 
cillating behavior of Bessel functions J n (7) 7 the hopping 
(2) 

rates n n can become complex- valued and the Hamilto- 
nian H2 , correspondingly, ceases to be non-Hermitian. 
This is clearly shown in Fig. 2 (a), where [tZn^] 2 becomes 
negative at a few lattice sites near n = 0. This circum- 
stance indicates that, for such indices, the hopping rates 
are purely imaginary. In spite of the non-Hermiticity of 
H2, its spectrum remains real- valued and BOs with the 
same period Tb as that of the WS Hamiltonian are found 
[20|. Fortunately, a second application of the intertwin- 



ing operator technique, assuming as a second energy level 
\±2 of Ho one of the two WS resonances adjacent to /ii, 
i.e. \±2 — ±F [corresponding to the quantum number 
I = ±1 in Eq.(19)], the Hermiticity of Hz is restored. 
As an example, Fig. 3 shows the behaviors of [ft^] 2 and 

= Vn^ — corresponding to [X2 — —F and for 

(3) 

two values of 7. Note that, for n — >• ±00, n n — » 1 and 
AVn^ settles down to constant values [AV^ — > for 
n — » +00, AVn^ — > 2F for n —> —00], so that the synthe- 
sized lattice described by H3 behaves like the WS (defect- 
free) lattice for site indices far from n = 0, 1, i.e. far from 
the lattice sites of WS localized states removed by the 
intertwining operator technique. As F is decreased, i.e. 
the localization length of the WS states is increased [see 
(20)], the localization length of the defect in the lattice is 
increased, as one can see by comparing Figs. 3 (a), (b) with 
Figs. 3(c), (d). The persistence of BOs in such synthesized 
lattices, corresponding to either single-site excitation at 
t = or to a multiple-site excitation with a broad Gaus- 
sian wave packet, is demonstrated in Figs.3(e-h). In the 
figures, the numerically-computed evolutions of \c n (t)\ 2 
versus t, as obtained by solving the Schrodinger equa- 
tion for the Hamiltonian H3 with the initial condition 
c n (0) = <5 n ,o (f° r the breathing BO modes, Figs. 3(e) and 
(g)] and c n (0) = exp[— (n - 10) 2 /64] (for the oscillatory 
BO modes, Figs. 3(f) and (h)] are depicted for two values 
of F. 

Tight-binding lattices with engineered hopping rates and 
energy sites corresponding to the ones shown in Fig. 3 
could be realized using arrays of evanescently-coupled 
optical waveguides with engineered size and distance, in 
which the distances between adjacent waveguides control 
the hopping rates n n whereas the channel widths (or re- 
fractive index changes) of the guides set the values of the 
site energies E n (see, for instance, [3, [H, Hl| and refer- 
ences therein). 

The technique of intertwining operators can be iterated 
by removing additional resonances from the WS ladder 
of Ho- Extended numerical simulations show that the 
resulting Hamiltonians turn out to be Hermitian pro- 
vided that couples of adjacent WS resonances are re- 
moved, whereas in the other cases the hopping rates can 
become imaginary at some lattice sites. As a general 
rule, the removal of a new WS resonance couple at ener- 
gies E = IF, (l±l)F introduces in the lattice new defects 
localized at around the lattice site n = I. For the Hamil- 
tonian H2S+1 obtained by removing from Ho s couples 
of adjacent WS resonances, one has ft! 2s+1 ^ — >• 1 for n — >• 
±00, AV; (2s+1) -> for n -> +00, and AF n (2s+1) 2sF 
for n — >• — 00. As an example, Fig. 4 shows the behav- 
ior of hopping amplitudes and site energy offsets for a 
synthesized lattice obtained by removal of the four WS 
resonances with energies E = 0, —F, — 25F, — 26F. As the 
number of removed WS resonances (and hence of defects 
in the lattice) increases, the behaviors of hopping ampli- 
tudes n n and site energies V n become highly irregular, 
as shown in the example of Fig. 5. Here the Hamiltonian 



5 



His is synthesized by successive application of the inter- 
twining operator technique that removes from the WS 
ladder spectrum the 12 resonances E = 42F, 41F, 26F, 
25F, 10F, 0, -F, -11F, -12F, -30F, -3LF. It is 
remarkable that, in such a rather irregular lattice with ex- 
tended defects, BOs still persists and exact self-imaging 
is attained, as shown in Figs. 5(c) and (d). 
As a final comment, it should be noted that the hop- 
ping rates and site energies of the synthesized lattices 
described by the Hermitian Hamiltonians H5, H7, 
... depend on 7, i.e. on the amplitude of the forcing 
F entering in Ho- Therefore, a change of the forcing 
parameter F in the original WS Hamiltonian Ho gives 
different lattice realizations for H3, H$, H7, ... (see, for 
instance, Fig. 3). This means that, as an additional term 
Hp = —f n \ n ) to Ho simply changes the BO period (just 
because the external forcing is changed from F to F + 
the addition of H p — —fn\n) to Hs-, H5, H7, ... destroys 
the onset of BOs. This is shown in Fig. 6 (a), where the 
temporal evolution of the site occupation probabilities 
|c n (t)| 2 is shown for the Hamiltonian His of Fig. 5, with 
an added perturbation term H p = —fn\n) for F = 0.2 
and / = —0.1. For comparison, the evolution of site oc- 
cupation probabilities for the lattice with the WS Hamil- 
tonian Ho perturbed with H p = —fn\n) (for the same 
values of F and /) is depicted in Fig. 6(b). Therefore, 
while in a WS ladder Hamiltonian BOs are observed for 
any value of forcing F (a change of F corresponds to a 
change of the BO period), in the synthesized lattices H3, 
H5, H7,... a change of the forcing is detrimental for the 
onset of BOs. Yet, it is remarkable that at a fixed forc- 
ing strength BOs can be observed in tight-binding lattices 
with defects, and even in greatly irregular lattices (like 
the one shown in Fig. 5). 



IV. CONCLUSIONS 

In conclusion, in this work it has been shown that sin- 
gle and multiple defects in a tight-binding Wannier- Stark 
lattice can be introduced such that BOs and the self- 
imaging property of the WS lattice are not destroyed. 
Such lattices are synthesized from the defect-free lat- 
tice by the technique of intertwining operators gener- 
ally employed in supersymmetric quantum mechanics to 
engineer the spectrum of Hermitian Hamiltonians. The 
energy spectrum of the synthesized lattices differs from 
the Wannier- Stark ladder of the defect-free lattice be- 
cause of the missing of pairs of resonances in the lad- 
der, thus ensuring the persistence of BOs. It is envis- 
aged that the lattice engineering technique proposed in 
this work could be extended to other coherent dynamical 
regimes, such as dynamic localization in presence of an 
ac (time-periodic) force [22| . It is also envisaged that the 
possibility to realize non-Hermitian tight-binding lattices 
with real- valued energies -mentioned in this work- would 
deserve further investigation and could stimulate novel 
studies in the framework of the rapidly developing field 
of non-Hermitian quantum mechanics [20|, l23| . 
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FIG. 3: (color online) (a-d) Behavior of the square of hopping 
amplitudes [ftL 3 ^] 2 and of site energy offset AVn 3 ^ = Vn^ + Fn 
for the Hamiltonian H3 synthesized from the WS Hamiltonian 
Tio by removal of the WS resonances fii = and \i2 — —F 
(see the inset at the top of the figure). Parameter values are 
k = 1 and F — 0.6 in (a) and (b), and k = 1 and F — 0.2 
in (c) and (d). (e) and (f): persistent BOs (breathing and 
oscillatory modes) in the the synthesized lattice corresponding 
to F — 0.6. (g) and (h): Same as (e) and (f), but for the 
synthesized lattice with F — 0.2. 
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FIG. 4: (color online) Same as Fig. 2, but for the Hamiltonian 
Tib synthesized from Tio by removing the four WS resonances 
E — 0,—Fj — 25F, — 26F. Parameter values are k — 1 and 
F = 0.6. 
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FIG. 5: (color online) (a) and (b) Same as Fig4., but for the 
Hamiltonian Wis synthesized from Ho by removing twelve 
WS resonances (parameter values are k = 1 and F = 0.2). In 
(c) and (d), the persistence of breathing and oscillatory BO 
modes, corresponding to either single-site excitation c n (0) = 
8n,o [(c)] and to a broad Gaussian wave packet excitation 
c n (0) = exp[— (n — 10) 2 /64] [(d)], are demonstrated. 
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FIG. 6: (color online) (a) Destruction of BOs (breathing 
modes) for the lattice Hamiltonian I-L12, of Fig. 5 due to the 
addition of the perturbation 7i p = —fn\n). Parameter values 
are k = 1, F = 0.2 and / = —0.1. In (b) the BOs with varied 
period, corresponding to the addition of 7i P to the WS lattice 
Hamiltonian Tio, are shown for comparison. 



